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CHAPTER 12 
THE REAL-TIME IMPLEMENTATION 


12.1 Introduction 

In this chapter, the final form of the real-time implementation of 
the PARC algorithm in a full duplex mode will be presented. The whole 
implementation will be decomposed into four parts. Each part is 
described in detail in the following sections. Block diagrams, flow 
charts and timing diagrams are given to help reader understand the 
implementation. The program listings are in Appendix G. 

Before presenting the implementation, the system components used 
will be described. A MAP-300 (Macro Arithmetic Processor) produced by 
Computer Signal Processing, Inc. is used to implement the speech coding 
algorithm (see Table 12.1). A host computer and MAP-300 system block 
diagram is shown in Fig. 12.1. The MAP-300 consists of six programmable 


processing elements as well as memory and peripherals: 


A. A Central Processing Unit (CSPU) functions as the executive con- 
troller by interpreting commands from the host computer, setting 
up and scheduling the other processors, and performing arithmetic 
and logic operations. 

B. An Arithmetic Processor (AP) consists of two subprocessors. 

An Arithmetic Processing Unit (APU) carries out the floating 
point arithmetic calculations. An Arithmetic Processing Scroll 
(APS) is a data addressing device for the APU. 

C. An Input/Output Scroll (I10S-2) transfers bit streams into or out 

of a modem. 


D. An Analog Data Acquisition Module (ADAM) samples and quantizes 


input analog signals into digital signals and stores them into 


main memory. 


ij Qty. Number 


12 


13 


Table 12.1 
MAP-300 Speech Processing Svstem 


Model 
1 1030 
1 2030 
1 2050 
1 2203 
1 2120 
1 3100 
1 4020 
2 4040 
1 5120 
1 5130 
1 6100 
1 6200 


1 03-1360585 


so Description = 


MAP-300 Processor 

8K x 32 MOS Memory 

16K x 32 MOS Memory 

8K x 32 MOS Memory 

2K x 32 Bipolar Memory 
PDP-11 Interface 

Model 2SM/ 1/0 Scroll 
Bus Switch 

Analog Data Acquisition Module 
Analog Output Module 
Expansion Chassis 
Auxiliary Power Supply 


GTE Speech Processing Interface 
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E. An Analog Output Module (AOM) converts input digital signals 
into analog signals and outputs them to a telephone module. 
F. A Host Interface Scroll (HIS) acts as a communication and data 


transfer medium between the host computer and the MAP-300 svstem. 


Tnis array processor system provides two classes of software. The 
first class is the software designed specifically for array processing. 
It is an executive-driver-subroutine package resident in the MAP-300 
and in the host computer which permits direct calling of routines to per- 
form arithmetic operation, to manage MAP memory, to define even higher 
level routines. A MAP user can execute array processing operations with 
simple Fortran calls in the host computer. The sophisticated programmer 
even can implement his own processing algorithms in MAP assembly 
languages. To implement the PARC digitization system, some new array 
and non-array functions, which are described later, are employed. The 
second class consists of utility programs including a cross-assembler 
and a cross-simulator. These may be used to add new routines to the 
array processing library. 

The array processor can be used as a peripheral to a host computer, 
using host peripherals for data storage and interface for data and 
commands transfer between them. However, after initialization, it can 
stand alone by using proper control functions in the host. 

The underlying design principles of the real-time implementation 
are as follows: First, in order to study system performance, the flexi- 
bility of changing system parameters such as speech algorithm parameters, 
buffer sizes and buffer starting locations must be provided. Second, in 


order to be in the real-time mode, the arithmetic execution time and the 


d overhead time have to be reduced. The arithmetic execution time can be 


I 20 
i 
decreased by efficiently utilizing the arithmetic processors. The over~ 
i ead time can be decreased by utilizing all processors as asynchronously 
and as simultaneously as possible. However, the necessary synchroniza- 
i tion between processors has to be carefully considered. 
| After several months of study, the final real-time PARC communica- 
tion system consists of the following programs: 
| 1. Host programs: Fortran PARC main program 
Fortran PARC host support program 
A 2. AP programs: PARC-transmitter pitch computation program 
| PARC-transmitter speech digitization program 
PARC-receiver 
d 3. CSPU programs: Encoder program 
Decoder program 
| Synchronization program 
| Transmitter buffer pointers update program 
Receiver buffer pointers update program 
| 4. I0S programs: ADAM program 
ADM program 
| I0S-2 program 
[ The main function of the host Fortran programs is to initialize the 
MAP-300 system with the real-time PARC algorithm, to set up the proper 
| command sequences and to initiate them. After that, the MAP-300 system 
can execute the PARC without any further commands from the host computer. 
| The two PARC-transmitter programs and the PARC-receiver program require 
| sophisticated arithmetic operations and will be executed in the AP. 
Those programs requiring logic and bit operations include encoder, de- 
5 coder, synchronizer and address pointers update programs. These will 
i 
a ey ee a ee ge a 
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be executed in the CSPU. A/D converting, D/A converting and bit-transfer 
will be executed in the ADAM, the AOM and the I0S~2 respectively. These 
five processors operate in parallel. The host Fortran modules will be 
described below; they present the whole picture of the PARC. The rest of 
the programs will be presented in the following sections. 

The flow chart of the real-time PARC main program is shown in the 
Fig. 12.2. The module can be divided into two steps: An initialization 
step and a processing step. In each step, each processor has its own 
role. 

1. The initialization step: 

In this step, the host computer will initialize the whole algorithm 
by reading in system parameters, configuring logical buffers, initializ- 
ing logical buffers, loading ADAM and AOM, and creating appropriate com- 
mand sequences called function lists. The logical buffers and the func- 
tion lists, which are fairly complicated, will be described later. In 
this step, only three processors in the MAP-300 system have been used as 
follows: 

A. The CSPU acts as the resource controller. Responding to a host 
command, it loads an A/D program to the ADAM, loads a D/A program 
to the AOM, configures the logical buffers, and loads appropriate 
arithmetic functions into the AP. 

B. The AP executes these arithmetic functions which will reset the 
logical buffers. 

C. The HIS acts as an interface between the host computer and the 


MAP-300. 


The Fortran host support program, which is the host support module 


for those new array and non-array functions in the real-time PARC | 


READ SYSTEM 
PARAMETERS 
CONFIGURE 

LOGLCAL BUFFERS 


INITIALIZE INITIALIZATION 
LOGICAL BUFFERS STEP 
LOAD 
ADAM & AOM 
CREATE 
FUNCTION LISTS 
PROCESSING 
STEP 


Fig. 12.2 The Flow Chart of the Main Program 
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algorithm, interfaces the main program to the host/MAP driver module. 
The main functions of this support program are as follows: 
1) Check each argument for the proper range and return. 
error code. 
2) Pack arguments for the MAP in a Function Control Block (FCB) 
3) Determine any host absolute addresses that must be evaluated 
while constructing the Data Control Block (DCB). 


4) Pass DCB to the host/MAP driver. 


2. The processing step: 

After initialization, the MAP-300 can execute the real time PARC 
algorithm without any further commands from the host computer. Because 
of the fixed delay for the whole svstem (refer to Section 2.2), a fixed 
time slot has been chosen as a main parameter to synchronize all the 
processors in the MAP-300 system as shown in Fig. 12.3 and Fig. 12.4. 

In each time slot, a particular function list, which is described below, 
will be executed. This fixed time slot is set to be 19687.5 microseconds 
which is precisely the period for ADAM to process 126 samples, the 
period for AOM to process 252 upsampled data, and the period for I0S-2 
to process 189 bits. 

The function lists, whose relationships are shown in Fig. 12.5, are 
eight sets of command sequences. When the real-time PARC is initiated, 
the first function list, which starts the ADAM, the AOM and the I0S-2, 
is executed. Afterward, the main control function list (the second 


function list) is executed by the CSPU. This function list assures the 


synchronization of the receiver. The execution of the list is as follows: 
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FUNCTION LIST 3 
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1. It initializes the receiver SHAT buffer and the AOM buffer. 


Also, it resets the receiver address pointers. 

2. If the channel synchronization does not exist, it keeps 
executing the synchronization operation until a channel 
synchronization is set. The synchronization operation 
consists of the PARC-transmitter, the encoder and the 
synchronizer. 

3. If the channel synchronization exists, it keeps executing the 
normal PARC operation until a channel synchronization is lost. 
The normal PARC operation consists of the PARC-transmitter, 


the encoder, the decoder and the PARC-receiver. 


4. Go tol. 

Because of parallel data processing which speeds up the whole system 
(refer to the Section 2.2), the double buffering technique is employed. 
Therefore, the synchronization operation and the normal PARC operation 
each consists of three function lists: A control function list, a 
function list which operates on the first half of the double buffers 
and a function list which operates on the second half of the double 
buffers. The control function list checks a semaphore which resides in 
the MAP memory and then executes a function list which will operate on 
the proper buffers. Th- reason for this is to have a continuity of data 
flow for the PARC-transmitter. For instance, suppose the receiver loses 
the channel synchronization after the normal PARC operation on the first 
half of the double buffers. Then the main control function list executes 
the synchronization operation. The control function of the synchroniza- 
tion operation checks the semaphore and then executes the synchronizatio 


operation on the second half of the double buffers. So, the transmitter 
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gets a continuous data flow. The main control function list will keep 
executing until the MAP-300 system is terminated. 

In the synchronization operation, which consists of the function 
lists 3, 4 and 5 (shown in Fig. 12.3), the role of each processor in a 
time slot is described as follows: 

1) The ADAM samples and quantizes the input analog signal from 
the telephone module and stores these samples into MAP main 
memory. In a time slot, it produces 126 quantized samples. 

2) The AOM acts as a D/A converter and outputs 252 silence 
samples since no information flows into the receiver. 

3) The IO0S-2 acts as an input/output interface between the 
MAP-300 and the channel modem. In a time slo-, it 
transfers 189 bits out of the encoder bit buffer and 
reads 189 bits into the decoder bit buffer. 

4) The CSPU functions as follows in sequence: 

A: It loads the AP with the pitch computation program. 

B: It executes the encoder. However, upon request from 
an APU interrupt, it suspends the current operation 
and 

C: Loads the AP with the speech digitization program. 
It restarts the encoder after the loading. 

D. It executes the channel synchronization program. 

E. It updates the PARC-transmitter address pointers. 

5) The AP executes the arithmetic part of the PARC-transmitter. 
After finishing each program, it will interrupt the CSPU to 
indicate that it is free for the next operation, it any. In 


this operation, the AP has two jobs: 
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A: The pitch computation program. 

B: The speech digitization program. 
These five processors function in parallel. The ADAM, AOM and ICS-2 
operate continuously without any interrupt. However, the CSPU and the 
AP have to finish their own job inside the time slot in order to be in 
real-time. 

The normal PARC operation (shown in Fig. 12.4) which consists of 
the function lists 6, 7 and 8, has the same ADAM and 10S-2 operation as 
the synchronization operation. However, the CSPU, the AOM and the AP 
have different operations as follows: 

1) The ADAM samples and quantizes the input analog signal 
from telephone module and stores these samples into MAP 
main memory. In a time slot, it produces 126 quantized 
samples. 

2) The AOM acts as a D/A converter and outputs 252 unsampled 
reconstructed speech samples to the telephone module. 

3) The I0S~2 acts as an input/output interface between the 
MAP-300 and the channel modem. In a time slot, it transfers 
189 bits out of the encoder bit buffer and reads 189 bits 
into the decoder bit buffer. 

4) The CSPU functions as follows in sequence: 
A: It loads the AP with the pitch computation program. 
B: It updates the receiver address pointers. 
C: It executes the encoder. However, upon request from an 

APU interrupt, it suspends the current operation and 


D: Loads the AP with the speech digitization program. It 


restarts the encoder after the loading. 


E. It executes the decoder. Upon request from an APU 
interrupt, it suspends the current operation and 
F. Loads the AP with the PARC-receiver program. It 
restarts the decoder after the loading. 
G. It updates the transmitter address pointers, 
5) The AP executes the PARC-transmitter and the PARC-receiver. 
After finishing each program, it will interrupt the CSPU 
to indicate that it is free for the next processing, if any. 
In this operation, the AP has three jobs: 
A: The pitch computation program. 
B: The speech digitization program. 
C: The PARC-receiver program. 
In this operation, the timing control is very important which will be 
discussed in the next section. 

The subsequent sections will describe each of the subsystems in the 
real-time PARC system. The next section explains the design of the 
PARC-transmitter. The complicated buffer system, buffer control and the 
speech digitization will also be described. The corresponding noiseless 
source encoder will be explained in Section 12.4. 

Section 12.3describes the design of the PARC-receiver. Although 
some parts of the system are the same as in the transmitter, the differ- 
ent design of the buffering and upsampling will be depicted. The 


corresponding noiseless source decoder will be explained in Section 3.5. 
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12,2 The PARC-Transmitter 

In this section, the main line of the real-time PARC-Transmitter, 
which resides in the AP, will be presented. The corresponding encoder 
will be discussed in Section 12.4. The block diagram of the PARC 
algorithm is shown in Fig. 2.2 and discussed in Chapter 2. 

Because of the limit on the size of the AP memory, the whole PARC- 
Transmitter algorithm cannot be implemented as a single AP program. So, 
the PARC-Transmitter has to be cut into two subprograms, namelv, the 
pitch computation program and the speech digitization program, as shown 
in Fig. 12.6. Those portions of the algorithm which operate on blocks of 
speech samples will be put in the pitch computation program. Those 
portions include: input gain factor calculation, system noise reducer, 
adaptive filter and pitch extractor. These elements of the algorithm 
which process speech samp le-by-sample will be combined into the speech 
digitization program. These elements are the adaptive quantizer, the 
inverse quantizer, the adaptive predictor and the pitch extraction. 

Before any further description of the design of this PARC- 
Transmitter, the buffer system has to be depicted. There are seven 
buffers employed by the PARC-Transmitter as shown in Table 12.2. 

Among them, ADAM, LEVEL and BIT buffers are double buffers which 
allow one processor to fill one half of the buffer while another 
processor processes on the other half. The SAMPLE, VHAT and SHAT 
buffers are circular buffers which allow continuous access. However, 
the additional address pointers, which indicate the starting locations 
of those circular buffers at the beginning of a time slot, have to be | 
considered. The PARAMETER buffer is a single buffer which stores of 


system parameters such as predictor coefficients, quantizer output 
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scaling factors, quantizer expansion factors, and the number of bits 
associated to each quantizer level. 

The PaRC-Transmitter employs four processors: The ADAM, the AP, 
the CSPU and the I0S-2. In order to achieve parallel processing, the 


double buffering technique is used. The function of these buffers 


een ed 


associated to the processors at a time slot is shown in Fig. 12.7. Ata 
time slot, the following processes are executed in parallel: 
1. While the ADAM is filling a half of the ADAM buffer, the AP 
is accessing the other half. 
2. While the AP is filling a half of LEVEL buffer, the CSPU is 
accessing the other half. 
3. While the CSPU is filling a half of the BIT buffer, the 
I0S-2 is accessing the other half. 
In the next two subsections, the pitch computation program and the 
speech digitization program will be described. The complicated buffer 


controller will also be presented. 


12.2.1 The Pitch Computation Program 


The pitch computation program contains the gain controller, a 
noise reducer, a buffer controller, the adaptive filter and the pitch 
y extractor. The APU flow chart and the corresponding APS flow chart, 
including the controlling flows are depicted in Figs. 12.8 & 12.9. The detail 
functions are described as follows and the parallel processing is shown: 
1. In response to the function list, the CSPU loads the 
APU and the APS modules of the pitch computation 
| program into the APU and the APS respectively. 


2. The CSPU then sets flag RI which will initiate the 


, | (| APS module. 
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The APS module sets flag RA which will start the APU module. 
The APS produces the addresses of the noise reducing parameter, 
the gain factor, 126 input data addresses from the ADAM buffer 
and 126 output data addresses to the SAMPLE buffer. 

The APU reads 126 input samples whose addresses are given by 
the APS, reduces their noise, multiplies by the gain factor 

and outputs to the SAMPLE buffer. The input samples from the 
ADAM are 11-bit precision plus a sign bit in a range of 
[2047/2048,-2048/2048]. After being read by the AP, these 


are converted into 32-bit floating point numbers. 


The APS checks the fullness of the SAMPLE buffer, sets the 


flag G2 if it is over the pitch repetition threshold. 


The APS again checks the fullness of the SAMPLE buffer. If it 
is under the filtering threshold, the APS sets the flag GO 

and jumps to Step 7. If it is over the threshold, the APS sets 
the flag AF1 and produces addresses for the filtering. The 
APU waits for the signals from the APS. If the flag GO is set, 
it goes to Step 7. Otherwise, it executes the filtering 


process which i3 described in Section 2.3. 


sane 
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10. 


The APS produces addresses for the data used to compute the 
pitch period. The APU computes the pitch period using the 
AMDF technique (refer to Section 2.4). In order to speed up 
the whole real-time system, all processors have to be utilized 
as asynchronously and as simultaneously as possible., i.e., 

the minimum use of the communication flags between the APU and 
the APS. However, because of the limit on the number of the 
registers in the APS, three flags are used in this portion to 
search for the pitch period as shown in Fig. 3.9. In this 
block, the APU determines how many iterations are left to be 
executed. The APS has to wait for the information, the flags 
AFO and AFl, at the end of each iteration. 

The APS produces addresses for the data used to calculate the 
pitch correlation coefficient. If the flag G2, which indicates 
the pitched repetition, is set, it produces addresses starting 
a pitched repetition block behind. The APU calculates the pitch 
correlation coefficient. 

The APS produces output addresses. The APU encodes the pitch 
period and the pitch correlation coefficient and outputs them 


to the MAP memory. 


The AP interrupts the CSPU to indicate that it is finished. 


APS 


Addresses 
for Pitch 
Period 

Searching 


. APU 
| 
| | Search for 


_ Clear AF3 
Update T 


Set AF2 |~ ~~~ — a 


} Fig. 12.9 The Flow Chart of the Pitch Period Searching Loop 


288 


12.2.2 The Speech Digitization Program 


The speech digitization program contains the pitched repetition, 
the adaptive predictor, a maximum number of processed samples con- 
troller, a reciprocal function, the pitch extraction, the adaptive 
quantizer, the inverse adaptive quantizer, the bit buffer controller, 
the underflow controller, and a parameter update. The APU and the APS 
flow chart including the control flows are depicted in the Fig. 12.10 
because the number of samples which are processed by the transmitter 
varies in each time slot, the flags AFO, AF3 and Gl are used to com- 
municate between the APU and the APS. The flag APO, which is controlled 
by the APU, is used to indicate the beginning of the digitization loop. 
The flag Gl, which is controlled by the APU and the APS, is used to 
indicate an underflow of the sample biffer or reaching the maximum rumber 
of samples permitted in the time slot. The flag AF3, which is controlled 
by the APU, is used to indicate that at least 157 information bits have 
been generated. 

The detail functions are as follows: 

1. In response to the function list, the CSPU loads the APU and 
the APS modules of the speech digitization program into the 
APU and the APS respectively. 
2. The CSPU then sets the flag RI which will initiate the 
APS module. 
3. The APS module sets the flag RA which will start the APU module. 
4. The APS produces the addresses of the system counters. The 
APU reads in the system counters. The system counters are the 


BIT counter, SAMPLE counter and the LEVEL counter. Because 


the fixed time slot is used as the main parameter for the | 
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real-time system, the I0S-2 transfers 189 bits in a time slot. 
However, as described in Chapter 2, only 157 bits are available 

to encode the quantizer levels. the BIT counter is used to 
control the PARC-Transmitter so as to generate just more than 

157 information bits, The SAMPLE counter is used to count the 
total number of samples processed by the PARC-Transmitter during 

a time slot. The address pointer update program will use this 
information to update the address pointers in the transmitter 
circular buffers. The LEVEL counter is used to limit the maximum 
number of samples which can be processed in a time slot. This 
assures real-time operation. 

The APU and the APS check the flag G2 which is controlled by the 
pitch computation program to indicate the condition of the pitched 
repetition. If flag G2 is clear, there is no pitched repetition 
and the APU and the APS go to Step 6. Otherwise, the APS produces 
the addresses of SHAT buffer for the pitched repetition. The APU 
executes the pitched repetition. The size of the pitched repeti- 
tion is a system parameter which can be input from the main Fortran 
program. The design of the real-time system assumes that the over- 
flow of the sample buffer can be absolutely controlled by the 
pitched repetition. Thus, the size of the repetition has to be 
carefully considered. 

The APS produces addresses for the predictor. The APU executes 
the predictor which employs a fourth order prediction. 

The APU checks the LEVEL counter and sets the flag Gl, if the 
counter is over the maximum number of samples allowed to be 


executed by the PARC-Transmitter. 
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10. 


12. 


13. 


The APU executes the reciprocal function to compute 1/ (RMS)? and 
1/o. To do the reciprocal of a 32-bit floating point number, 
the APU usually employs an 8 x 256 read only memory which has 
only 8-bit precision in the mantissa part. Thus, a special 
subalgorithm is used to increase the precision. 

The APS produces addresses for the pitch extraction. The APU 

executes the pitch extraction. 

The APS checks the SAMPLE buffer and sets the flat Gl if the 

buffer is empty. 

The APU and the APS quantize the input sample. Three flags are 

employed to communicate between the APU and the APS. The way 

to quantize an input sample is as follows: 

A. The APU signals the APS to produce the threshold address. 

B. The APU compares the magnitude of the input sample with the 
threshold. 

C. If the threshold is larger than the magnitude, the APU 
signals the APS to stop the searching and jumps to the 
inverse quantizer portion. Otherwise, the APU checks 
whether the threshold is the largest one. If it is, the 
APU signals the APS and jumps to inverse quantizer 
portion. If it is not, control is returned to Step A. 

This design allows variable quantizer levels and reduces searching 

time, because most of the input samples are dropped in the first 

two levels. 

The APU and the APS execute the inverse quantizer. 


The APU updates and checks the BIT counter and sets the flag 


AF3 if the 157-bit is reached. 


14. The APU and the APS update the predictor parameters. 
i 15. If the flag Gl, which indicates an underflow of the SAMPLE 
buffer or reaching the maximum number of samples permitted 
in the time slots, is set, the APU goes to Step 16. If the 
flag AF3, which indicates the condition of BIT counter, is 
set, the APU goes to Step 17. Otherwise, the APU sets the 
flag AFO and goes to Step 6. The APS waits for the signals 


from the APU. If the flag AFO is set, it goes to Step 6. 


| 

| 

| Otherwise, it goes to Step 17. 

| 16. If this step is executed, there has been an underflow in the 

SAMPLE buffer, i.e., less than 157 bits are generated. The 

| APU updates the BIT counter to account for a NULL code which 
indicates underflow. If the BIT counter is still less than 
157 bits, the BIT counter is reset to synchronize the counters 
between the PARC-Transmitter and the encoder. 

17. The APU and the APS output the contents of the system counters 


and parameters. 


18. The AP interrupts the CSPU to indicate that it is finished. 


12.3. The PARC Receiver 
| In this section, the PARC~Receiver, which resides in the AP, will 
be presented. The corresponding decoder and channel synchronizer will 
be explained in Section 12.5. The principal elements of the receiver, 
which is shown in Fig. 2.8, are also part of the transmitter. So, the 
implementation of receiver is much the same as that of the transmitter. 
Before the description of the design of the PARC-Receiver, the 


buffer system has to be depicted. There are six buffers employed by 


, is the receiver as shown in Table 12.3. 
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Among them, the AOM and the LEVEL buffers are double buffers which 
allow one processor to fill one-half of the buffer while another proces- 
sor to process on the other half. This increases the parallel ability 
of the receiver. The VHAT, the SHAT and the BIT buffers are circular 
buffers which allow continuous access. However, the additional address 
pointers, which indicate the starting locations of these circular buf- 
fers at the beginning of each time slot, have to be considered. The 
structures of the LEVEL buffer and the BIT buffer will be explained in 
Section 12.5. The PARAMETER buffer is a single buffer which stores svstem 
parameters such as predictor coefficients, quantizer output scaling fac- 
tors and quantizer expansion factors. 

The flow chart of the receiver is shown in Fig. 12.11. The detailed 
functions are described as follows and the parallel processing is shown: 

1. In response to the function list, the CSPU loads the APU and 
the APS modules of the PARC-Receiver program into the APU and 
the APS respectively. 

2. The CSPU sets the flag RI which will initiate the APS module. 

3. The APS then sets the flag RA which will start the APU module. 

4. The APS produces 126 input addresses of the SHAT buffer and 

252 output addresses of the AOM buffer. The APU processes the 

upsampling operation on 126 input reconstructed speech samples 

and produces 252 nnsampled data points. The upsampling operation 

employs the linear interpolation technique which is explained in 

Section 2.8. The APU also limits the values of output data to 

be in the range of [2047/2048,-2048/2048] which is acceptable 


for the AOM. 


ee 


a ae = Er 


298 


APU APS cSPU 


Load 


APU & APS 


‘-@-> 


Addresses * 
Upsampling for ©) 
Process Upsampl ing i 
Process v 


Read Systen Addresses . 
Parameters for System 6) 

and Parameters 
Counter & Counter 


Pitched 
Repetition 


Fig. 12.11 The Flow Chart of the PARC-Receiver Program ed 


ge ert 


mime He eB eRe eel ommelCUeDlUCUCUlUlC GUC Ul eel eel eel eee 


Set AFO 


Fig. 12.11 


APU 


Read First 
Quantizing 
Level 


Adaptive 
Predictor 


Inverse 
Adaptive 
Quantizer 


Enpitched 
Process 


Update 
Predictor 


aYra>»n 


yes 


Output 
System 

Parameters 
& Counter 


I 
| 


Set AF3 |- 


--AFO,_AF1 _ 


eecer-- HH 


The Flow Chart of the PARC-Receiver Program 


APS 


Address of 
First Data 


Addresses 
for 

Adaptive 

Predictor 


Addresses 
for 

Inverse 

Quantizer 


Addresses 
for 
Enpitched 

Process 


Addresses 
for 
Predictor 


Parameters 


ves 


Addresses 
of System 
Parameters 
& Counter 


<Q) 


aS 


<- 


¢-©- 2 


emt rept aae. 


300 


10. 


11. 


The APS produces addresses of the system parameters and a 
counter. The APU reads the system parameters and the counter. 
There is only one system counter, the SAMPLE counter, which 

is used to count the total number of samples processed by the 
PARC-Receiver during a time slot. The receiver address pointer 
update program will use this information to update the address 
pointers of the receiver circular buffers. 

The APU checks the first data in the LEVEL buffer and decides 
the condition of the pitched repetition. If no pitched 
repetition is employed, it sets the flag AFl and goes to 

Step 7. Othrrwise, it sets flag G3 and processes pitched 
repetition. The APS waits for signals from the APU. If the 
flag AFl is set, it goes to Step 7. If the flag G3 is set, 

it produces addresses for the pitched repetition operation. 

The APS produces the second data which is the first received 
quantizing level in the LEVEL buffer. The APU reads the first 
quantizing level. 

The APS produces the addresses for the adaptive predictor. The 
APU processes the adaptive predictor which employs a fourth 
order prediction. 

The AP operates the inverse adaptive quantizer. Two flags, AFO 
and AFl, are used to communicate between the APU and the APS. 
The APS produces the address of V(K-T). The APU restores the 
pitch redundancy as shown in Eq. (2.1). 


The APS produces addresses of the predictor parameters. 


The APU updates the predictor coefficients. 
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12. The APU reads the next data in the LEVEL buffer and checks 
H whether it has reached the end of the buffer. If it has, 
the APU sets the flag AF3 and goes to Step 13. Otherwise, 
i it sets the flag AFO and gces to Step 8. The APS waits for 
signals from the APU. If the flag AF3 is set, it goes to 
Step 13. If the flag AFO is set, it goes to Step 8 to 


continue the receiver loop. 


eS eed 


13. The AP wutputs the system param ters and the counter to the 


eee 


i MAP memory. 
14. The AP interrupts the CSPU to indicate that it is finished 


with the PARC-Receiver. 
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12.4 Noiseless Source Coder 


This section describes the implementation of the noiseless source 
coder on the central processor of the MAP 300 called CSPU. As discussed 
in Section 2.7, the source coder performs several functions. Its primary 
purpose is encoding the quantizer levels and the side information for each 
block of Ne samples processed by the PARC transmitter. In addition, it 
supplies tue bit pattern for frame synchronization and the parity code 
for channel error control. 

At its input, the encoder shares a double buffer containing the pitched 


repetition indicator and N, quantizer levels each with the PARC transmitter. 


B 
The corresponding pitch reduction parameters 8,T are in two sets of locations 
in the scalar table. At its output, it shares a double buffer containing 
189 bits each with the I0S, the processor which interfaces with a digital 
modem. 

The encoder strives to form a frame of 189 bits from the quantizer 
levels and the side information generated by the transmitter. The frame 
is subdivided into three Hamming blocks of 63 bits each. A Hamming block 
contains 57 information bits and 6 parity bits. Each of the 57 information 
bits in a block is protected against errors. Depending on whether or not 
pitched repetition is indicated, the output frame assumes one of two formats 
Figs. 12.12 and 12.13. The encoder proceeds sequentially to build the different 
fields in the output frame. It handles one Hamming block at a time. The 
parity codeword for the block is updated edch time a bit is output. Immediately 


after the 57th information bit in a Hamming block is transferred, the parity { 


codeword for that block is output. The codeword is then initialized for the 


next block, and the encoder continues to output further information bits. This 
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B Pitch period T 

c Pitch correlation coefficient 8 
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Fig. 12.12 Normal Format for Bit Frame Output of Encoder 


omar = 


Hamming Biock 1 


Field Contents 
A Synchronization bit 
B Pitch period T 
c' False 8, for pitched repetition 
Cc Pitch correlation coefficient 8 
D1 Quantizer levels 
El Parity bits for error control 


Fig. 12.13 Format of first Hamming Block during pitched repetition 
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process is repeated for three Hamming blocks, completing the output 
frame of 189 bits. 

Fig. 12.14 shows the flow chart for the encoder. The program listing 
is contained in Appendix G. The encoder starts with field A, Fig. 12.12. 
It fetches the previous sync bit from a location where it was saved and in- 
verts it. This new value is saved for use in the next frame. It is also 
output to field A. Next the encoded value of T is fetched and transferred 
to the output buffer, field B. The value of the pitched repetition indicator 
is then checked. If affirmative, the code 0000000 is transferred to field 
C*, Fig. 12.13. Nexe the partially encoded value of 8 is fetched. It must 
be between 0 and 96 indicating one of 97 possible levels. Using this as an 
index, a 7-bit code is retrieved from the 8 encode table and output to field 
C, Figs. 12.12, 12.13. After this the quantizer levels are fetched one by 
one. The corresponding code-lengths and code-words are retrieved from the 
Q-level encode table and transferred to the output buffer until the fields 
D1, D2, and D3 are filled. This should coincide with the exhausting of the 
Np Q-levels generated by the transmitter. 

An exception to this occurs when the sample buffer runs out. The block 
of quantizer levels do not generate enough bits to fill up the fields Dl, 
D2, and D3. In this case, a null code 11111110 is output to indicate the 
end of quantizer level information. If there is still more space available, 
it is padded in with 1's. The null code and the padding bits are error- 
protected the same as any other information bits. 

The source code used here is a variable-length to variable-length mapping, 
and it is not necessary that fields Dl, D2, and D3 are exactly filled by 


encoding an integer number of quantizer levels. The extra bits after a 
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Hamming block is full are passed to the appropriate field in the next 


Hamming block. 


All the programs in each processor need to be able to process a block of 


information within a time frame of 19687.5 yp sec for the algorithm to 


keep in real time. This necessitated giving program efficiency priority 


over flexibility. The following is a list of the constraints on the encoder 


program as a result of this: 


(i) All buffers used by the encoder must be located in Bus 1. They 


can be relocatable, as long as they do not occupy memory locations above 64K. 


(ii) All the inputs to the encoder program must be in the fixed number 


format. The quantizer levels must be pre-multiplied by 2 to facilitate 
direct indexing. The values of 8 and T must be partially encoded. T must 
be limited to 64 values from 0 to 63, and 8 to 97 values from 0 to 96. 
Pitched repetition and the end of quantizer levels must be indicated by 
negative numbers. 

(iii) The source code has the following constraints. The code for 
level 1 is 0. The locations for the run length code and the null code ia 
the Q-level encoding table are 0 and 24, respectively. The other 11 codes 
for the quantizer levels Q are at locations 2Q. 

(iv) Because the LSB of a code-word is transmitted first, the code- 
word received at the decoder is inverted, e.g., 110 + 011. To compensate 
for this, the code-words stored in the Q-encoding table must be inverted 
and right justified. 

(v) The Q-encoding table must contain (L-1) for the various code 
lengths, L. 


(vi) The code to indicate pitched repetition is 0000000. The 97 


}-- 
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quantized values of 8 may not use any pattern with 2 or less 1's in it. 
This is to give the pitched repetition indicator additional error protection. 
(vii) The run length value can be selected at initialization. 
(viii) Finally, the block of Np quantizer levels should not generate 
more than 196 bits. Since the encoder does not monitor the maximum number 


of bits generated, the transmitter must ensure this limit. 


. 


Serene 


12.5 Synchronizer, Decoder 

This section describes the implementation of the decoder on the central 
processor of the MAP 300 called CSPU. The decoder essentially inverts at 
the receiver the operation performed by the encoder at the transmitter. 

From each received frame of 189 bits, it generates a block of Ve quantizer 
levels, the pitched repetition indicator, and the pitch extraction parameters 
for use by the PARC receiver. In addition, it performs several other operations. 
It monitors transmission errors, and corrects up to one error per Hamming 
block in the received frames. Because the received bit stream is arranged 

in frames, the decoder acquires the initial frame synchronization using the 
synchronization bit pattern transmitted by the encoder. Later, it monitors 
the received sync pattern to ascertain that synchronization is not lost. 

It monitors the state of the sample buffer in the PARC receiver to make sure 
the buffer does not overflow or underflow. Finally, it controls the mode of 
operation of the entire PARC algorithm. 

The different operations outlined here are performed in different 
modules in the decoder. The decoder program is subdivided into three modules: 
the initialization module; the synchronization acquisition module; and the 
decoder module. The function to be performed and correspondingly the module 
to be used is determined by the three values of a function select switch, 

S!. The first time the program goes in for sync acquisition, several parameters 
and buffers need to be initialized. This mode is indicated by a value of 0 

for Sl. Subsequent sync acquisition operations are indicated by a value of -l 
for Sl. After frame sync has been established, the decoder can perform its 

task of inverting the bit stream to quantizer levels. This mode is indicated 


by a value of 1 for $1. 
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The decoder shares a double buffer with the PARC receiver at its out- 
put. A second function select switch, $2, indicates which half of the 
double buffer the decoder must use in the current time frame. 

At its input, the decoder shares a circular buffer 1024 bits long with 
the IOS. The decoder expects to receive 189 bits in each time frame. 
However the clock that controls the inflow of the bit stream is at the 
other end of the communication system. It is different from the clock that 
measures the time frame at the decoder. On the average, the decoder does 
receive 189 bits per time frame. However, because of the different clocks, 
there could be temporary deviations. The large circular buffer which is 
several times the frame size allows for these deviations without causing a 
conflict between the decoder and the IOS over the bits being read out of 
and written into the buffer. The IOS, which receives the bit stream from 
the digital channel, has its base pointer which indicates the position of the 
base of the current frame to be entered into the buffer. When the PARC 
algorithm is initiated, the base pointer of the decoder which indicates the 
base of the frame of bits it should process in the current time frame is set 
several blocks behind the input pointer controlled by the I0S. Each process- 
or updates its pointer independently at the end of the time frames. As long 
as the input and output rates match on the average, no conflict should occur 
and the two pointers should remain reasonably separated. 


On being called, the decoder checks S1 and transfers control to the 


ee 


appropriate module. The module performs its function, and then it sets up 


‘ A 
oer} 


the control and updates the relevant parameters and information for the 
operation in the next time frame. The program listing of the decoder is { 


included in Appendix G. The following subsections describe in detail the 


operation of the three modules that make up the decoder. i 
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12.5.1 Synchronization Acquisition Module 

Before the decoder can start inverting the received bit stream, the 
boundaries of the frames, marked by the sync bit, must be located. This 
task is achieved by the synchronization module, over several time frames. 
The sync acquisition algorithm is described in Chapter 3. The implementation 
here varies slightly from the description in Chapter 3. This is to reduce 
the amount of computation and memory storage required and to make the pro- 
gram more time efficient. The block diagram of the sync acquisition module 
is shown in Fig. 12.15. 

The encoder inserts an oscillating bit Ss; in field A, Fig. 12.12, for 
synchronization. The decoder generates a similar pattern of expected sync 
bit values, S.- During sync acquisition, the correlation of S and each of 
the 189 bit positions starting from the base pointer is computed. A cumulative 
correlation tally from one time frame to the next for each of the bit positions 
is stored separately. Just in case the sync sequence being generated at the 
decoder is the inverse of the sync pattern being received from the transmitter, 
a cumulative tally of correlations with a is also saved for each of the bit 
positions. The bit position corresponding to the maximum cumulative tally 
in each time frame is kept. When the correlation tally for a particular bit 
is maximum for ten consecutive frames, it is assumed the sync bit has been 
located. That position, marking the beginnings of frames generated by the 
encoder, is saved. 

If at the end of a sync acquisition operation, it is determined that 
sync has not yet been acquired, the program sets up for one more sync 
operation in the next time frame. The base pointer in the circular buffer 
is updated by 189. The current value of S, is inverted and saved. And the 


i 
function select switch S2 is changed to indicate opposite buffer affiliations 
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for the next time frame. 


If synchronization is achieved, preparations are made for subsequent 


Orman 4 


decoder operations. In addition to the above controls and parameters, several 


Se mmaaial 


others are updated. Function select switch S1 is changed to 1 to indicate 

a decoder operation in the next frame. Although no decoder operation has so far 
been performed, a dummy buffer consisting of 126 level 1's is prepared for 

use by the PARC receiver. Error check is conducted on the first Hamming 


block of the next bit frame; and up to one error is corrected. T, 8 and the 


output later. This leads to the location of the first bit in the field Bl, 
Fig. 3.12, of the next frame. This is called the lst information bit, 


| pitched repetition indicator for the next frame are decoded and saved for 
i information about quantizer levels. Its location is saved for use in the 


next frame. The cost function for sync monitor is initialized to -1. And 


aA 


i the expected sync value S is set to match the next received sync value 


itl 
| Siete Finally, the state of the receiver sample buffer is initialized to 


a buffer-full condition. The control is then returned to the MAP executive. 


12.5.2 Decoder Module 
This module is the heart of the decoder. It corrects for transmission 
errors, if possible; monitors frame synchronization; inverts the received 
bit sequence to information usable by the PARC receiver; and performs buffer 


control on the receiver sample buffer. Its block diagram is shown in Fig.12.16. 


i 
| 


After it decides its buffer affiliation by checking $2, it does error 
detection and correction on the 2nd and 3rd Hamming blocks of the current bit 
frame and the lst Hamming block of the next bit frame. The reason for this 
offset between the bit frame and the error control frame is to make the 
extra bits generated during the encoding of the current frame available to 


the decoder when it performs the decoding. The extra bits reside at the 
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beginning of field Bl, Fig. 12.12 of the next frame. 


The next step is sync monitor. The basis of the algorithm is explained 
in Chapter 3. The sync variable starts with an initial value of -1, FFFF in 
hexadecimal representation. The received sync bit Sy is compared with its ex- 
pected value S, using the logical operation EXOR. If equal, the bits in the 
sync variable v, are shifted right one position. Bit 15, the sign bit is re- 
tained 1 in this operation. If not equal, the sync variable is shifted left two 
positions, causing 0's to be inserted in the least significant bits of V5° 

Bit 11 of v5 is monitored. If it is 0, sync is declared lost. The module 
sets the switch Sl to 0 to indicate a series of sync acquisition operations 
starting in the next time frame, and returns control to the MAP executive. 

As long as bit 11 of Vy is 1, sync is declared preserved, and the module 
continues to its next phase of operation, source decoding. The pitch reduction 
parameters 8 and T which were decoded and saved in the previous time frame are 
output to their appropriate locations in the scaler table. The pitched rep- 
etition indicator, also decoded in the previous time frame, is transferred to 
the output buffer. If affirmative, the receiver buffer pointer is updated by 
the length of the repetition block. 

After this, the pitch reduction parameters and the pitched repetition 
indicator for the next time frame are decoded and saved. This leads to the 
position of the lst information bit in the next frame. If there were any 
extra bits in the current frame, the decoder will know where to look for 
them when it needs them. 

The source decoding is a tree searching operation along the coding tree. 
Each successive bit received represents a node in the tree. Starting from 
the root of the tree, bits are read in, until the last bit received represents 


the termination of the tree. The corresponding quantizer level(s) are determined 
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from the decoding table and transferred to the output buffer. This process 
continues until either a null code is received or all the bits in the current 
frame are exhausted. At the end of the the decoding operation, the location 
of the lst information bit for the next frame is saved. 

The decoder can correct transmission errors to a point, after which 
errors filter through to the received bit stream. Errors have two effects 
on the decoding procedure. The quantizer levels generated are wrong. And 
the number of quantizer levels Ny in the current block can be different 
from what it should have been. This causes the state of the sample buffers 
in the transmitter and receiver to lose correspondence. To correct for 
this, the decoder provides a steady drift in the state of the receiver buffer. 
When N, for a received block is 126, it indicates that the transmitter was 
processing silence when it generated this block and its sample buffer was 
empty. Correspondingly, the receiver buffer should be full. A few extra 
level 1's are added to the q-buffer generated for the receiver; so that if 
errors have caused the receiver to move away from its buffer~-full condition, 
it should slowly drift back. This alleviates the mismatch between the 
buffer states at the receiver and the transmitter. 

Finally, the state of the receiver sample buffer is checked. If it 
underflows, extra level 1's are added to remedy it's negative state. If it 
overflows, the extra samples are deleted. 

After this, the module updates various controls and parameters for 
another decoder operation in the next time frame. It then returns control 


to the MAP executive. 


12.5.3 Initialization Module 


The initialization module precedes the first of a sequence of sync 
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acquisition operations. Sync acquisition requires two buffers initialized 
for its proper operation. This is done by the initialization module. The 
locations which contain the bit positions with max correlation tallies 
from the last ten frames are initialized to -1. The buffer where the 
correlation tallies are stored is initialized to 0. The switch S1 is set 
to -1 to indicate subsequent sync acquisition operations. 

In addition, the length of the pitched repetition block is retrieved 


from the function control block and saved for later use in the decoder. It 


‘also gets the base address of the input buffer. It then transfers control 


to the sync acquisition module. 


12.5.4 Decoder Constraints 

As with the other modules in the PARC algorithm, the stress in developing 
the decoder program is execution efficiency rather than flexibility. This 
is to try and meet the rather tight limitations of real time operation. The 
following constraints are imposed on the decoder: 

(i) All buffers used by the decoder must reside on Bus 1. They can 
be relocatable, as long as they do not occupy memory locations above 64K. 

(ii) The input for this program is a circular buffer of length 1024. 
The information bit must be contained in the LSB of each word. 

(iii) All the outputs of the decoder are in the fixed point number 
format. The 8 and T are partially encoded and are represented by 7 and 6 
bits respectively. 8 is allowed one of 97 codewords, 0 to 96, and T one of 
64, 0 to 63. Pitched repetition and end-of-quantizer-levels are indicated 
by negative numbers. 


(iv) The run length value can be selected at initialization. 
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12.6 Program Timing and Speed 


The programs in the real time implementation of PARC get executed in 
parallel on the different processors of the MAP 300. The set of programs 
on a given processor must be executed once in each time frame. The number 
of samples processed by these programs is not constant. It varies from 
about 50 during the voiced regions of speech to about 520 during the 
transition regions following voiced speech. For the algorithm to operate 
at its peak performance, the programs should be able to process the maximum 
number of 520 samples within a time frame. This is, however, limited by 
processor speeds and program efficiency. 

The average number of samples that need to be processed by any program 
is 126; and it is imperative that the programs be able to process at least 
a few extra samples if necessary. Currently, the programs do satisfy this 
minimum requirement. The pitched reduction program, PARC transmitter, and 
PARC receiver which operate on the AP can process about 150 samples per 
time frame. The encoder and decoder programs which operate on the CSPU can 
process about 200 samples per time frame. So currently, the limit on the 
number of samples processed is set at 150. This exceeds the minimum require- 
ment only marginally, and speeding up the programs would improve the performance 
of the algorithm considerably. Some ways to achieve this improvement are 
suggested below, and it is expected that different programs will execute 20% 
to 50% faster as a result of these modifications. 

The pitched reduction program which currently requires about 5 msec can 
be speeded up 100% by doubling the parallelism in its computation. The trans- 
mitter and receiver programs need to be reorganized in terms of their register 


usage and operation sequencing. This should provide approximately 30% speed 


banner al 


meant 


improvement. These changes should increase the number of samples processed 
on the AP by about 50%. 

The encoder can be speeded up 25% by rearranging the encoding table 
and changing all memory access from indirect to direct addressing. This 
would require the input quantizer levels to be 2(q-l) instead of 2q as thev 
are now. The decoder can be improved slightly by rearranging its decoding 
table and operation sequencing. 

Two of the three buses on the MAP 300 have slower MOS memories. Changing 
these to the fast bipolar memory would also help the speed on the various 
processors because of the large number of memory transf2rs performed in each 
time frame. With these modifications, it is expected the algorithm would 
be able to process 250-300 samples per time frame instead of the 150 it does 
now. Although this still allows the algorithm to operate well below its 


peak level of fidelity, this is about the level of performance that can be 


expected given the currént speed of the MAP 300 array processor. 


APPENDIX E 
RESAMPLING PROGRAM 


During the course of this project, because of the different sampling 
rates at which digitized speech was available from various sources and the 
different rates at which it was utilized, it was found necessary to develop 
an efficient resampling program to enable us to easily change sampling rates 
of speech. The following briefly describes the resampling algorithm and the 
operational details of the program. Fig. E.l1 shows the flow chart, and 
program listings are included at the end of this appendix. 

Sampled speech is the pulse amplitude representation of an analog speech 


signal at a sequence of time points separated by T At another sampling 


1° 
rate, it is the PAM representation of the same signal at a different set of 


points separated by T The new set can be obtained by interpolation from 


2° 
the first set of samples. So as also to limit the bandwidth of signals, 
interpolation was performed by discrete convolution of the original samples 
with the truncated impulse response of an ideal low pass filter. A first 


order Hamming window was used for truncation to reduce the effects of truncation. 


The impulse response and the Hamming window have the following form: 


H, (t-t) = 2F sinc (2F mr t) -» < T < @ 
and Wy (t-t) = a+ (l-a) cos (n t/2T) -T<1t<T 
a = 0.7 


The convolution function is the product of these two. The discrete 
version of the convolution function is scaled down by Ths to compensate for 


the scaling introduced by discrete filtering. 
C, (nT,.-(m-i)T, ) wet + (l-a) cos [((mT,-nT,) + iT,)/ NT.J])} 
i 2 1 T, 1 2 1 r 


{sinc[2Fn((mT-nT,) + iT,)]}, 


en ne ee ee 


—* 


1) Name of input file 
2) Name of output file 
3) Comments on output file 


4) Output Sampling rate, R> 


Open input file 
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1) Length of file, Ly 
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Fig. E.1 Flow Chart for the Resample Program 


Loop: 
till input 


file is exhausted 


oats 


Here nT, identifies the nth output sample and mT, the closest priveding 


1 
input sample. 


Because the phase (mT, - nT,) changes arbitrarily, the correlation 


i 
coefficients need to be computed for each output sample. However, if th: 

two sampling rates are integer submultiples of a higher frequency ! Ee tallied 
the reference frequency, the C5 need to be evaluated only once at the referees 
frequency. The phase (nT, - nT? is always a multiple of ex ind a subset 

of these coefficients can be utilized to perform the convelution. This re- 
duces the computation required in the program by about an order of magnitude. 
For the different rates required in this project, the frequency 96000 is in 
integer multiple and was selected as the reference frequency. 

Two versions of the resampling program were developed. Ia the first. 
all computation, initialization as wel] as the convolution (which comprises 
the bulk of the computation in resampling) is done on the PDP 11/60. In the 
second, the convolution is performed on the MAP 300 array processor using the 
SNAP-II software package of array functions supplied by CSPI. The PDP 1! 
does the initialization and handles the I/O interfacing between the disk and 


the MAP 300. The second version is almost an order of magnitude faster than 


the first. 


E.1 Operating Details 


The program is segmented into four modules. The main program acquire - 
all relevant parameters required for the resampling operation, either ‘ror: 
the header record of the input speech file or from the terminal. This 
program prompts for each input from the terminal, and it requested, vives 
a brief description, shows an acceptable example, and specifies the beund- 
on the value of the input. The second module performs the convoiution for 
resampling. The other two modules assist in the initialization bv computing 
the convolution coefficients. Between the two versions of the program only 
the second module is different. The appropriate module is selected during 
the process of task-building. 

To obtain an executable task image, the appropriate modules must be com- 
piled and task-built. The following describes this process for the RSXI1IM 
operating system on a PDP 11 computer. The compilation step for each module 


has the following form. 
>FOR MODULE = MODULE/14 


The task-building step is slightly different for the two versions. The 


command format for the first version is: 


>TKB 
TKB> QRESAM/CP/FP = MOD1, MOD2, MOD3, MOD4 
TKB> / 

ENTER OPTIONS: 

TKB> UNITS = 7 

TKB> ACTFIL = 7 

TKB> // 


The underlined text is the prompt by the computer. The following text 
is the response by the user. 
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For the second version, the SNAP library must also be included in 


in the task-building step. 


>TKB 
TKB> MRESAM/CP/MP = MOD1, MOD2, MOD3, MOD4, SNAPLIB/LB 
TKB> / 


ENTER OPTIONS: 
TKB> UNITS = 10 
TKB> ACTFIL = 10 
TKB> // 


The program requires the input speech file to conform to the Notre 
Dame speech file format. Fig. E.2 illustrates this format. The first record 
is the header. The seven fields in it are sentence identification number, 
sampling rate, number of samples in the files, lower and upper cutoff 
frequencies, truncation number for the convolution window, and comments 
identifying the speech file. The upper cutoff frequency is not used in this 
program. The truncation number which specifies the length of the convolution 
window is internally computed by the program based on the input sampling rate. 
Each of the subsequent records in the file contain 16 speech samples in the 
15 format. The amplitude of speech samples is limited between -2048 and 2047. 
The speech file must end with a null record. 

The output file generated by the resampling program also conforms to 


this format. 


First i i : H i i 
Record: Header A AZ} Az 'A4 ! AS 4 AG ' A7 

} ‘ ' ¢ t t 

' ' ' ' ‘ 


Subsequent 16 Speech Samples 
Records: Speech IO ci I i A EE NS | 


| 16 Speech Samples 


Last Null 
Record: End Eee at ee ee 


Format for each record of 80 columns 


Header: (615, 10X, 2042) 
Speech: (1615) 


Fig. E.2 Notre Dame Speech File Format 
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APPENDIX F CVS) ALGORITHM 


In this appendix, the CVSD algorithm will be discussed. Before the 
design of this real-time CVSD system can be presented, it is essential 


to understand the structure of CVSD systems. 


F.1 The CVSD System 


In this section, the structure of a CVSD system is presented. The 
system described here is identical to that of [1]. A schematic of an 
encoder and decoder pair for CVSD is shown in Figure F.1. 

A CVSD encoder is simply a delta modulator with an adaptive stepsize. 
In order to minimize the difference between adjacent samples, the input 
signal s(k) is normally oversampled. Typically, the sampling rate is 16 
k samples per second. Since one bit per sample is transmitted, the trans- 
mission rate is 16 K bits per second. 

The encoder quantizes the difference e(k) between an input speech 


] sample and its estimate which is predicted by a first order predictor. 
e(k) = s(k) - a * s(k-1) 


where s(k-1) is the reconstructed speech of the input speech sample at time 


k-1. Then the encoder outputs a single bit, b(k), for each corresponding 
input sample where 
if e(k)>0 


b(k) = 
if e(k)< 0 


The adaptive stepsize logic derives its output A(k) from the bit stream 


b(k) and an appropriate initial state, 


A(k) = B * A(k-1) + g(k) 


4upOoul (G) Uapoouy (¥) waysdg (SAD oUL Td ‘Sta 
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(A Can ear If b(k) = b(k-1) = b(k-2) 
g(k) i max 


A * (1 - 8) otherwise 
min 
For a long run in which no three consecutive bits are identical, 


the stepsize approaches Anin’ On the other hand, if the many consecutive 


bits are the same, the stepsize will approach Oe age 


The estimate of the difference signal is given by 
e(k) = b(k) * A(k) 
Then the reconstructed speech is computed from 
s(k) = a * s(k-1) + e(k) 


Because the reconstructed speech in the transmitter is determined 
solely by the bit stream b(k) and an agreed initial state, the receiver can 
produce the same reconstructed speech samples if no channel error occur. 

Normally, there is no best criteria to judge the performance of a 
speech algorithm. However, the basic signal to noise ratio criteria is 


used here. 


M 
) sty 
M 


SNR = 10 * 10819 est 


Y ( s(k) - 3(k) )? 
k=1 


F.2 The Real-Time CVSD System 


In this section, the design of the real-time CVSD system is presented. 
The major design decisions and real-time algorithm are discussed. The 


program listings are followed. 


F.2.1 Design Overview 


In order to design a real-time speech system, several points have to 


be considered. First, the flexibility of chanzing system parameters suchas 


ememet ——_ 25 2B EB GE GF 5 a ee | 


| 
l 
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speech algorithm parameters, buffer sizes and buffer starting locations 
must be provided. This allows system performance to be studied without 
any further modification of the real-time program. However, this 
approach will increase overhead time. Therefore, a "pre-bound" concept 
will be introduced. Second, the arithmetic execution time has to be 
minimized. This can be done by efficiently utilizing the arithmetic 
processors. Therefore, the "hiding" technique which hides additions 
inside the period of multiplication will be employed. Third, in order to 
decrease overhead time, all processors have to be utilized as asynchronously 
and as simultaneously as possible. However, the necessary synchronization 
between the processors has to be considered. 

The real-time CVSD implementation employs five MAP processors to do 
speech processing asynchronously and simultaneously. The whole algorithm 
is divided into two steps: An initialization step and a processing step. 
Each processor has its own role in each step. 

1. The Initialization Step 

In this step, the host computer will initialize the MAP-300 system, 
set up the appropriate command sequence, and initiate the real- 
time CVSD algorithm. Each processor in the MAP-300 has to do 

the following functions: 

A. The CSPU acts as the resource controller. Responding to a 
host command, it loads a data acquisition program into the 
ADAM, loads digital-to analog conversion program into the 
AOM , configures the buffers, loads an arithmetic function 
into the AP which will reset the containts of the buffers, 
builds a recurrent sequence of commands called a function 


list and prebinds the APS module of the CVSD program. 
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B. The AP executes the arithmetic function which is loaded by 
the CSPU. The purpose of this execution is to reset the 


buffers. 


2. The Processing Step 

After initialization, the MAP-300 can execute its functions 

without any further commands from host computer. The role of 

each processor are described as follows: 

A. The CSPU acts as the resource controller. Responding to 
each interrupt, it updates the containts of the correspond- 
ing processor control block. Responding to the function 
list, it checks the availability of the appropriate buffer 
and initiates AP operation by loading the CVSD program 
from MAP memory. The synchronization is also done by CSPU 
which is described in the next section. 

B. The APU executes the real-time CVSD algorithm. 

C. The APS produces data addresses and acts as the controller 
of the real-time CVSD system. 

D. The ADAM samples the input speech signal from telephone 
module and stores them into MAP memory. 

E. AOM converts input digital samples into an analog signal and 


output it to telephone module. 


In the next two subsections, the detailed description of these two steps 
will be presented. 
F.2.2 The Initialization Step 

The initialization step, which is controlled by host computer, includes 
reading system parameters, giving commands to the CSPU and building a 


function list. The algorithm is shown in Fig. F.2. 


Read System 
Parameters 


Configure 
Logic Buffer 


Reset Buffer Initilaization Step 


Load 
ADAM & AOM 


Create 


Function List 


Start 
CVSD System 


Processing Step 


Fig. F.2 Algorithm for Initialization Step 
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There are two points worth to note as follows: 


1. 


Because of parallel processing in the MAP-300, the synchronization 
between the AP and the AOM, and between the AP and the ADAM has 
to be considered. Memory is the common place which every 
processor has to access. Therefore, the synchronization is 

done by checking buffer availability as follows: 

A If the first ADAM buffer is busy, CSPU goes to wait. Other- 
wise, goes to Step B. 

B. CSPU initiates AP. AP starts to read input samples from the 
first ADAM buffer, execute CVSD function and output results 
into the first AOM buffer. _ 

C. If the second ADAM buffer is busy, CSPU goes to wait, 
otherwise, goes to Step D. 

D. CSPU initiates AP. AP does the same function as Step B except 
it reads samples from the second ADAM buffer and outputs 
euatea ts the second AOM buffer. 

E. Go to Step A. 

This structure is included in a function list. 

In order to increase the flexibility of changing system 

parameters, there are some substitutive attributes inside the 

APS module. In a normal operation, CSPU gets values for those 

attributes from a host command and then binds them with the APS 

module before loading it into the APS memory. In order to 
decrease this overhead time, a "pre-bound” has to be done prior 

to real-time execution. The prebinding process will create a 


new APS module from old APS module by permanently binding the 
attributes. At the execution time, all normal operations of 


binding and buffer protection will be bypassed. Thus, the over- 


head time is decreased. 
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F.2.3 The Processing Step 
In the processing step, each of three processors has a program. Among 

them, the data acquisition and data output programs, which are included in 
the Simple Notation for Array Processing library (SNAP library), are fairly 
simple. The actual CVSD algorithm is executed by the AP processor. 

The real-time CVSD program can be divided into two parts as follows: 
1. APU module of CVSD algorithm which is shown in Fig. F.3 will do 

the whole arithmetic operation described in the Section F.1. The 

input and output data are directed by APS module. The "hiding" 


technique is employed to reduced execution time. 


te 


The APS module of the CVSD algorithm which is shown in Fig. F.4 
will produce address sequences of input and output data for APU module. 
Also, APS module acts as the controller of the real-time CVSD system. 
This AP control protocol is shown in Fig. F.5 and described as follows: 
* CPU starts APS-input by setting RI 
* APS sets the program counter for write routine and starts APS~ 
output by setting RO. 
* APS starts APU by setting RA. 
* After producing addresses for input data, APS-input turns itself 
off by clearing RI. 
* When FI is cleared, APU leaves the loop. 
* After producing addresses for output data, APS-output turns itself 
off by clearing RO. 


* APU turns itself off by clearing RA. 


* CSPU is interrupted by both RAI (RA off) and DNI (AP done). 
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Fig. F.4 APS Flow Chart 
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F.3 Conclusion 

The execution time for the real-time CVSD algorithm ia approximately 
12 micro-seconds per sample. So, this implementation can be used with a 
sampling rate up to 83 KH in the real-time mode. A set of parameters which 


gives a very good quality of speech is as follows: 


‘max = 750 
Ain = 25 
a = 0.94 
B = 0.99 


A two seconds sample speech is executed by this algorithm, and a 11.87 db 


signal-to-noise ratio performance is obtained. 
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{1] R. W. Schafer et. al., "Tandem Interconnection of LPC and 
CVSD Digital Speech Coders" Final Report to Defense Communications 
Agency, DCA Contract No. DCA 100-76-C-0073, November, 1977. 
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APPENDIX G 


PROGRAM LISTINGS FOR REAL TIME IMPLEMENTATION 
OF PARC ON MAP-300 


This appendix contains the program listings of all the modules in the 
real time implementation of the PARC algorithm on the MAP-300 array processor. 
There are twenty modules in all, some in the PDP~11, and others in the various 
Processors of the MAP 300. Ten modules execute in the MAP-300. Nine of these 
have host support programs in the PDP-11. And finally, there is the main 
program in the PDP--11 which sets up and initializes the real time PARC operation 
on the MAP-300. 

The following is a list of these programs and the processors they execute 


on. 


1. PARC Mainline PDP~11 

2-10. Host Support Subroutines PDP-11 

11. Receiver Update Module CSPU /MAP-300 

12. Transmitter Update Module CSPU/MAP-300 

13. Initialization Module CSPU/MAP-300 

14. Pitch Computation Module AP /MAP-300 

15. PARC Transmitter Module AP /MAP-300 

16. PARC Receiver Module AP /MAP-300 ; 
17. Encoder Module CSPU/MAP-300 { 
18. _- 10S Simulator CSPU/MAP-300 . 
19. Decoder Module CSPU/MAP-300 


20. Digital 1/0 Module 10S /MAP-300 
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Operator Sequence to Run PARC at DCA 

The two MAP's at DCA have different logical device names, MP and 
MA. There is a different task images corresponding to each MAP, called 
PARP.TSK and PARA.TSK. The task names for these programs are NDP and 
NDA. 

To run the real time PARC algorithm, the following steps must be 
followed for each MAP. For the steps shown below, the devices are 
called MX. The underlined text is the computer response at the terminal. 


The rest of the text is the operator input. 


> LOA MX: (Load the map device) 
> INS MXLD/TASK = ... XLO (Install map loader program) 
> XLO (Load the map executive including PARC routines. 
The executive is called NDQL) 
INS PARX (Install the initialization routines for PARC) 
> NDX (Run the initialization program) 


MODE ?(1 - Internal Loopback, 2-Full Duplex) 
2 (Mode Selection) 


These steps for setting up and executing the PARC algorithm are 
combined into an indirect command file called PARC.CMD. If the indirect 
command file is used, the operator only needs to run the program (NDX) 


and select the mode. 
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